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Abstract 

We study numerical integration on the unit sphere 8 2 C M 3 using equal weight 
quadrature rules, where the weights are such that constant functions are integrated 
exactly. 

The quadrature points are constructed by lifting a (0, m, 2)-net given in the 
unit square [0, l] 2 to the sphere S 2 by means of an area preserving map. A similar 
approach has previously been suggested by Cui and Freeden [SIAM J. Sci. Comput. 
18 (1997), no. 2]. 

We prove three results. The first one is that the construction is (almost) optimal 
with respect to discrepancies based on spherical rectangles. Further we prove that 
the point set is asymptotically uniformly distributed on S 2 . And finally, we prove an 
upper bound on the spherical cap L2-discrepancy of order AT -1 / 2 (log A) 1 / 2 (where 
A denotes the number of points). This improves upon the bound on the spherical 
cap ^-discrepancy of the construction by Lubotzky, Phillips and Sarnak [Comm. 
Pure Appl. Math. 39 (1986)] by a factor of ydogA. 

Numerical results suggest that the (0, m, 2)-nets lifted to the sphere S 2 have 
spherical cap L2-discrepancy converging with the optimal order of A -3 / 4 . 



1 Introduction 

We consider the unit sphere S> 2 = {z = (zi,z 2 ,z 3 ) G IR 3 : ||z|| = \/ z\ + z\ + z\ = 1}. 
Let / : § 2 — > R be integrable. Then we estimate the integral J g2 / da 2 , where a 2 is the 
normalized Lebesgue surface area measure (that is J* g2 do"2 = 1), by a quasi-Monte Carlo 
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type rule 



1 N ~ l 



N 

k= 



where Z N = {z , . . . , z N _i\ C § 2 are the quadrature points on the sphere. Since the 
surface area measure is normalized to 1, it follows that we have Qn(J) = J§2 / dcr 2 for 
every constant function /. In the following we review known results to put the result of 
this paper into context. Although some results are known for spheres S d of dimension 
d > 2, we only state them for the sphere S 2 since this paper only deals with the 2-sphere. 

Using Stolarsky's invariance principle [ID] (also see [10] and [II]). it follows that the 
worst-case error for numerical integration in a certain reproducing kernel Hilbert space is 
given by the spherical cap L 2 -discrepancy of the quadrature points. To obtain quadrature 
points we use a transformation from [0, l] 2 to S> 2 which preserves area. Specifically, we 
use the transformation $ : [0, l] 2 — > S> 2 with 

$(2/1,2/2) = ^2 cos(27ryi) ^y 2 -y\,2 sin(27ryi) ^ y 2 - y\, 1 - 2y^j . 

The function $ maps axis-parallel rectangles in the unit square to zonal spherical rect- 
angles of equal area. It is natural then that the discrepancy on the sphere with respect 
to spherical rectangles is the same as the discrepancy with respect to rectangles in [0, l] 2 . 
Since point sets Zn for which a discrepancy based on ^-regular test sets (K fixed, see 
Sjogren [37J) converges to zero as N — > 00 are uniformly distributed over the sphere, we 
show that the digital nets lifted to the sphere via $ are also uniformly distributed. 

Furthermore, we study the spherical cap L 2 -discrepancy of point sets obtained in this 
way. We prove that the order of convergence for (0, m, 2)-nets lifted to the sphere via $ is 
(9(A r_1//2 (log iV) 1 / 2 ). This improves upon the bound on the spherical cap /^-discrepancy 
in [26] of 0(N~ 1 ^ 2 log N) by a factor of \/log N. On the other hand, the optimal order of 
convergence is 0(N~ 3 ^ 4 ). But numerical experiments do suggest that the (0,m, 2)-nets 
lifted to the sphere via $ achieve optimal convergence rate. We conjecture that this order 
of convergence is the correct one. 

Due to the difficulties of having a satisfactory notion of 'bounded variation' on § 2 
there is no Koksma-Hlawka inequality on § 2 per se. However, the concept of uniform 
distribution of a sequence of N -point configurations with respect to every function in a 
function space, say Sobolev spaces over S 2 , can be used as quality criterion. For example, 
Cui and Freeden [T3| introduced the concept of a generalized discrepancy on S> 2 which 
involves pseudodifferential operators. Using the operator D = (— 2A*) 1 / 2 (— A* + 1/4) 1 / 4 
(A* is the Beltrami operator on § 2 ) of order 3/2 they arrived at 



N 

iir £/(**)- //d. 

k=i 



N 

' v 2 



< v / 6D(K...,x w };D) 



n 3/2 



where / is from the Sobolev space % 3 / 2 (S> 2 ). In this notion, a sequence {Z^} of iV-point 
configurations is called D-equidistributed in % 3 / 2 (§ 2 ) if liniAr^oo D({zi, . . . , Zjy}; D) = 0. 
The generalized discrepancy associated with D can be easily computed by way of 

1 - 

4tt [D({z 1; . . . , z N }; D)] 2 = 1 - _ ^ log (1 + ||z, - 2fc || /2) 2 . (1) 

kl=l 



2 



Sloan and Womersley [38] showed that [D({zi, . . . ,zn}; D)] 2 has a natural interpretation 
as the worst-case error for the equally weighted quadrature rule Qn associated with the 
points Zi, . . . ,Zjy for functions / from the unit ball in "H 3 / 2 (§ 2 ) provided with the re- 
producing kernel K(y,z) = 2[1 — log(l + \\y — z\\/2)]. In [TT] this approach is followed 
further, yielding an even simpler notion of discrepancy (see Section [3J, namely 

4 1 N 

[D({ 2l , . . . , z N })] 2 = - - — \\*t ~ *k\\ (2) 

k,e=i 

also used in the setting of the Sobolev space K 3//2 (S 2 ) now provided with the reproducing 
kernel K(y,z) = (8/3) — \\y — z\\. Low-discrepancy configurations in the above contexts 
can be found by maximizing the respective double sums in Q and Q. This leads into 
the realm of the discrete (minimum) energy problem on the sphere where points are 
thought to interact according to a Riesz s-potential l/r s (s > 0) or logarithmic potential 
log(l/r) (s = 0) and r is the Euclidean distance in the ambient space. It is known that 
the minimizer Zn (N > 2) of the associated s-energy functionals form a sequence which 
is asymptotically uniformly distributed for each fixed s > 0. We refer the interested 
reader to [HI El El d [HI H21 [201 I2H E21 [H] . We remark that Sun and Chen [JT] 
employ 'spherical basis functions' (as a counter part to radial basis function on spheres) 
to investigate uniform distribution on spheres. They also show that the minimizers of the 
functionals induced by a spherical basis function are asymptotically uniformly distributed. 

Let C(z,t):={y 6 S> 2 : (y,z) < t} be a spherical cap and let C:={C(z,t) : z G 
S 2 , — 1 < t < 1} denote the set of all spherical caps. In [ID], Stolarsky established a 
beautiful invariance principle on the 2-sphere (and higher-dimensional spheres), 

i N r r 

— J2\\ze-z k \\+4[D 2 (Z N ;§ 2 ,C)} 2 = / \\y - z\\ da 2 (y) da 2 (z), (3) 

k,£=l ■' ■ ■' ■ 

connecting the sum of distances, the /^-discrepancy (with respect to spherical caps 

COM)), 

1/2 



D 2 (Z N ;S 2 ,C):-- 



-i 



\z N nc(z,t)\ ( ... 

a 2 (C(z,t)) 



N 



da 2 (z) dt 



and the distance integral. We observe that, on S> 2 , the discrepancy in ^ is essentially (up 
to a factor 2) the L 2 -discrepancy. Originally, Stolarsky used his invariance principle and 
sharp result for discrepancy estimates by Schmidt [36J to establish bounds for the sum of 
distances. Beck [1] then improved Stolarsky's lower bound, finally arriving at 

N 

cN-V 2 < I I \\x-y\\da 2 (x)da 2 (y)-^J2\\z e -z k \\<CN- 3 / 2 (4) 

J k,i=l 

for some universal positive constants c and C independent of N. Consequently, rela- 
tions (|1J) yield lower and upper bound for the /^-discrepancy by means of the invariance 
principle which are sharp with respect to order of TV: 

c'jV-3/4 < d 2 (Z* n ;S 2 ,C) < C'N- 3/4 

for a sequence of optimal L 2 -discrepancy iV-point configurations Z* N on S> 2 . Observe, 
that optimal L2-discrepancy and optimal spherical cap discrepancy configurations satisfy 
estimates with the same order iV~ 3 / 4 (apart from a possible y/\og N term), see PJ. 
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1.1 Quasi-Monte Carlo rules in the unit square 

Quasi-Monte Carlo algorithms /(/) = 4 ^2^=o f( x n) are used to approximate integrals 
1(f) = L -a /(a;) da;. The crux of the method is to choose the quadrature points 
x , . . . , Xjsr-i as uniformly distributed as possible. The difference to Monte Carlo is the 
method by which the sample points Xq, . . . ,Xn-i £ [0, l) 2 are chosen. The aim of QMC 
is to chose those points such that the integration error 

N-1 



I f(x)dx - — J2f(x n 



achieves the (almost) optimal rate of convergence as iV — > oo for certain classes of func- 
tions / : [0, l] 2 — > R. For instance, for the family of functions / with bounded variation 
in the sense of Hardy and Krause (for which we write ||/||hk < °°) it is known that the 
best rate of convergence for the worst case error is 



e = sup 

/,II/IIhk<oo 



1 N ^ 

/ f(x)dx - — J2f{x r 



N- 1+£ for all e > 0. 



4M 2 n=0 

More precisely, there are constants c, C > such that 

cN^y/logN <e< CN^y/logN, 



see [To] . 

There is an explicit construction of the sample points Xq, . . . ,x^-i for which the 
optimal rate of convergence is achieved. One criterion for how uniformly a set of points 
Pn = {xq, ■ ■ ■ , Xn-i} is distributed in the unit square is the star discrepancy 



N-1 



D*(P N ;[0,1] 2 ,K*) = sup \5 PN (y)\ 

ye[o,i] 2 



y) = jjYl 1 "*6[o.») ~ Area ([°> v))> 



n=0 



where [0, y) = UUA^Vi) with 1/ = (Vi>V*) e [°> ^ n * = il ^) ■ V G t ' !] 2 }> 
Area([0, y)) = yiT/ 2 is the area of [0, y) and 



L a:;e[0,y) 



1 if a^e [0,y), 
otherwise. 



The quantity Sp N (y) is called the local discrepancy (of P/v)- 

The connection between this criterion and the integration error is given by the Koksma- 
Hlawka inequality 



f 1 N ~' 

/ f(x) dx - — V /(aj n 



<27*(P^r; [0, l] 2 ,K*)\\f\\ RK . 



Informally, a sequence of points x ,Xi,... G [0, l) s is called a low-discrepancy se- 
quence, if 



D*({x , . . . , a^Tv— i } j [0, l] 2 , ft*) = 0(N- L (\ogNY) as iV ^ oo. 
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Notice that for such a point sequence, the Koksma-Hlawka inequality implies the optimal 
rate of convergence of the integration error (apart from the power of the logiV factor), 
since for a given integrand / the variation ||/||hk does not depend on P N and N at all. 

The concept of digital nets introduced by Niederreiter [30] provides the to date most 
efficient method to explicitly construct points sets Pn = {xo, . . . ,xn-i} G [0, l) 2 with 
small discrepancy, that is 

D*(P N ;[0,1] 2 ,TV) ^CN- 1 log N. 
They are introduced in the next section. 

1.2 Nets and sequences in the unit square 

In this section we give a brief overview of (digital) (0, m, 2)-nets and (digital) (0, 2)- 
sequences. For a comprehensive introduction see [15] . 

The aim is to construct a point set Pn = {xo, . . . ,xn~i} such that the star discrep- 
ancy satisfies D*(P N ; [0, 1} 2 ,1Z*) < CN" 1 logiV. To do so, we discretize the problem by 
choosing the point set Pn such that the local discrepancy 5p N (y) = for certain y e [0, l] 2 
(those y in turn are chosen such that the star discrepancy of Pn is small, as we explain 
below). 

It turns out that, when one chooses a base b > 2 and N = b m , then for every natural 
number m there exist point sets P^m = {x , . . . , x b ™-i} such that 5p bm (y) = for all 
V = (2/1,1/2) of the form 

y i = ai /b dl fori = 1,2, 

where < dj < b dl is an integer and d\ + d 2 < m with d\, d 2 > 0. A point set Pn 
which satisfies this property is called a (0,m,2)-net in base b. An equivalent description 
of (0, m, 2)-nets in base b is given in the following definition. 

Definition 1. Let b > 2 and m > 1 be integers. A point set P\,™ C [0, l) 2 consisting 
of b m points is called a (0,m,2)-net in base b, if for all nonnegative integers d\,d 2 with 
d\ + d 2 = m, the elementary interval 

2 

n 

i=l 

contains exactly 1 point of Pym for all integers < a« < b di . 

It is also possible to construct nested (0, m, 2)-nets, thereby obtaining an infinite 
sequence of points. 

Definition 2. Let b > 2 be an integer. A sequence Xq,Xi, . . . 6 [0, l) 2 is called a (0,2)- 
sequence in base b, if for allm > and for all k > 0, the point set Xkb m , x^b^+i, • • • , X( k+ i) b m_ 1 
is a (0,m,2)-net in base b. 

It can be shown that a (0, m, 2)-net in base b satisfies 

771 

D*(P N ;[0,l} 2 ,n*)<C b — , 



Oi Oi + 1 \ 

Jjdi ' Q di J 
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and the first N points x , . . . , xn_i of a (0, 2)-sequence in base b satisfy 



D*({x ,...,x N ^};[0,l] 2 ,TZ*) <Q 



N 



for all N > 1, 



where C& > depends on 6 but not on m and iV. See [T5], [3T] for details. 

Explicit constructions of (0, m, 2)-nets and (0, 2)-sequences can be obtained using 
the digital construction scheme. Such point sets are then called digital nets (or digi- 
tal (0, m, 2)-nets if the point set is a (0, m, 2)-net) or digital sequences (or digital (0, 2)- 
sequence if the sequence is a (0, 2)-sequence). 

To describe the digital construction scheme, let b be a prime number and let Z& be 
the finite field of order b (a prime power and the finite field F& could be used as well). 
Let C\, C*2 € Z™ xm be 2 matrices of size m x m with elements in Z&. The ith coordinate 
x Uj i of the nth point x n = (^ n ,i,^n,2) of the digital net is obtained in the following way. 
For < n < b m let n = n Q + + • • • + n m _ 1 6 m_1 be the base b representation of n. Let 



n 



(n Q , 



n, 



m—li 



£ ZP denote the vector of digits of n. Then let 



For y nA = (y n> 



) T G Z™ we set 

__ l/n,?.,l 



Vn,; 



To construct digital sequences, the generating matrices G\,Ci are of size oo x oo. 

The search for (0, m, 2)-nets and (0, 2)-sequences has now been reduced to finding 
suitable matrices C\ y G-i- Explicit constructions of such matrices are available and where 
introduced (in chronological order) by Sobol [39], Faure [18], Niederreiter [30], Niederreiter 
and Xing j32j |33l H5], as well as others. For instance, to obtain a digital (0, m, 2)-net over 
Z2 one can choose 



Ci 



/ 1 


. 


. 


o\ 






(0 





1 ' 






















and 


c 2 = 









. 1 














. 


. 


1 J 






V 



(!) 



... (V) 

fm-2\ fm—V\ 

\m-2J \m-2J 

C:D J 



where the binomial coefficients are taken modulo 2. Notice that the matrices C\,Ci can 
be extended to 00 x 00 matrices yielding generating matrices for a digital (0, 2)-sequence 
over Z 2 . This example was first provided by Sobol' in [39] . 



2 Spherical rectangle discrepancy of point sets on the 
unit sphere S 2 

We introduce a transformation to lift (0, m, 2)-nets to the sphere S 2 . To do so, we represent 
each point on the sphere S 2 by using scaled spherical coordinates. We define T : [0, 1) x 
[0, 1] ^ § 2 by 

T(9, (f>) = (cos(27r6 l ) sin(7T0), sin(27r#) sin(7T0), cos(7r</>)) . 
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Let < 0i < < 2 < 1 and < 0i < 2 < 1. The part of the sphere 

= 0) G § 2 : X < < 2 , 0! < < 2 } 



has area 
Let 



Area(tt eiid2:(f>U( p 2 ) = 2vr(0 2 - 0i) [cos(tt0i) - cos(tt0 2 ) 



Area(Q eii e 2 ^ lj02 ) cos(tt0i) - cos(tt0 2 ) 



Area(§ 2 ) 

be the area of ^e 1 ,e 2 ,<t>i,<i>2 normalized with respect to the surface area of the unit sphere. 
Therefore we have r ,i,o,i = 1- We can now define a discrepancy measure of point sets on 
the sphere with respect to sets ^6>i,0 2 ,</>i,</>2- 

Definition 3. Let Z N = {z , . . . , z N _{\ C § 2 . Then the extreme spherical rectangle 
discrepancy of Z N on S> 2 with respect to Q = {^9 1 ,e 2 ,(f> 1 ,(t> 2 : < 0i < 2 < 1, < 0i < 2 < 
1} is defined by 



where 



D N (Z N ; S> 2 , Q) = sup 

O<6>i<0 2 <1 

O<01<</>2<1 



N-l 



iv, 



e 2 Ai,4>2 Oi,02,4>i,<h 



n=0 



1 i/z n ? ^i,8 2 A,fe 
otherwise. 



The spherical rectangle star- discrepancy of Z N on S> 2 with respect to Q* = {Qo,e,o,4> '■ 
< < 1, < < 1} is defined by 



D* N (Z N ;S\Vl*) = sup 

0<6K1 

O<0<1 

1 \r^N—l 



1 ^ 



N 



ZnGf!o,e,o,. 



0,6,0,, 



n=0 



The expression ^ J2 n =o lzetie^,*^ denotes the proportion of points of Z N in ftfc.fc.fc.te- 
Note that the set £V,0 2 ,</>i,<fe always includes the North Pole (0, 0, 1) and never the South 
Pole (0,0,-1). 

Further, the discrepancy measure includes the discrepancy of spherical caps centered 
at the North Pole (0,0,1) and South Pole (0,0,-1). The spherical cap at the North 
Pole is obtained by using Q ,i,o,<t>- Let > be small. Then -L Y.n^ l z n m ,i fiA is tne 
proportion of points in the spherical cap ^0,1,0,0 an d 47rr o ,i ) o,0 is the area of the spherical 
cap. Hence 



N-l 



Zn£f^0,l,0,. 



r o,i,o, 9 



n=0 



measures the discrepancy at the North Pole. The South Pole is included in the following 
way. Let < 1 be close to 1. Then § 2 \ f2 , 1,0,0 is the spherical cap centered at the South 
Pole. Then 

1 



N-l 

1-iVl 

N ^ 



n=0 



2 n Gf2o,l,0,. 



n=0 
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is the proportion of points in the spherical cap centered at the South Pole. Further, 
47r(l — r ,i 5 o,^) is the area of the spherical cap § 2 \ O Oj i,o,0- Hence 



Y ^b6^o,i,o,# — ^0,1,0,. 



N 



n=0 



N-l 



Zne§ 2 \f2o,l,0,</> 



Area(§ 2 \ Q . 



1,0,0; 



n=0 



Area(§ 2 



measures the discrepancy at the South Pole. 

In the following we construct point sets Zn = {z , . . . ,zn-i} on the sphere § 2 such 
that Dn(Zn] S 2 , Q) is small. This can be done by relating the spherical rectangle discrep- 
ancy to an analogous discrepancy on [0, l] 2 . 



Definition 4. Let P/v = {^cb • • • > a^jv-i} Q [0, l) 2 . Then the extreme discrepancy of P, 



n 



on [0, l) 2 with respect to 71 = {[ai,a 2 ) x [ci,c 2 ) : < ax < a 2 < 1, < Ci < c 2 < 1} is 
defined by 



where 



D{P N ;[0,iy,n) = sup 

0<ai<a 2 <l 
0<ci<c 2 <l 



L £c n e[ai,a2)x[ci,c 2 ) 



1 N ^ 

~jy 1 a=ne[ai,a2)x[ci,c 2 ) ~ ( a 2 ~ Oi)(c 2 - Ci) 



n=0 



1 i/cc n G [ai,a 2 ) x [c l} c 2 ) 
otherwise. 



The star- discrepancy of Pn on [0, l) 2 mtt respect to TV = {[0, a) x [0, c) : < a < 
l,0<c<l}zs defined by 



D*(P N ;[0,iy,K*)= sup 



0<a<l 
0<c<l 



iV-l 



1 ^ 

1 ^e[o,a)x[o,c) - ac 



n=0 



Let Pat be an arbitrary point set in [0, l) 2 . It is known, see for instance [HH Chapter 3], 



that 



D*{P N ; [0,l) 2 ,ft*) < D N (P N ; [0,l) 2 ,ft) < 4D*{P N ; [0,1) 2 ,^ 



The following lemma establishes a connection between the discrepancies on [0, l) 2 and 
§ 2 defined above. 

Lemma 1. Let P/v = {xo, • • • ,xn-i} £ (0, l) 2 and x n = (x ni i,x nj2 ) for < n < N . Let 

arccos(l — 2x n 2 ) 



z n T I x n i 



:2 



for < n < N 



and Zn = {zo, . . . , zn-i}- Then z n e S /or < n < N and 

D N {Z N - § 2 , Q) = D N (P N ; [0, l) 2 , K) 

and 

D* N (Z N ;S 2 ,n*) = D N (P N ; [0,1) 2 ,^*) 
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Proof. It is trivial to check that \\z n \\ = 1 for all n. To show the second part, we set 
Q\ = d, 9 2 = d2, 0i — [arccos(l — 2ci)]/V and 02 — [arccos(l — 2c 2 )]/tt. Then we have 

r #i ,0 2 ,0i,0 2 = («2 - — 2 ^ — ^ = («2 - «i)(c 2 - Ci). (5) 

Note that the points (0, 0, 1) and (0, 0, —1) are excluded from the set Zn since we assume 
that Pn C (0, l) 2 . The mapping from P/v to Zn is therefore bijective. Therefore we have 

l#„efi fll ,62,0!, 02 = l(y n ,l,(axCCOs(l— 2j/„ i2 ))/7r)6[01,02)x [01,02) = l(j/n,l.?/n,2)e[ai,a2)x[ci,C2)- 

Thus the remaining results follow from Definitions [3] and |4j □ 
Remark 1. ./Vote that sin(arccos(l — 2x)) = \Jx — x 2 , whence z n = (z n i, z n2 , z nj3 ) where 



z n>1 = 2 cos(27rx n ,i) a / x n>2 - x 2 n 



.,2' 



z n>2 = 2sin(27rx n> i)yx n) 2 - x 2 i2 , 

/or < n < N . To shorten the notation we define 

arccos(l — 2x 2 ) 



$(xi,x 2 ) = Tlxi 



71 



^2 cos(27rxi) \J x 2 — x 2 , 2 sin(27rxi) \J x 2 — x 2 ,, 1 — 2x 2 ^ 



for (xi,x 2 ) G [0, l] 2 . Hence, z n = §(x n ) for < n < b m . 
Further, for J C [0, l] 2 we set 

$(J) = {$(x) G § 2 :xeJ}. 



2.1 Digital nets on the sphere 

The previous lemma can now be used to obtain a bound on the spherical rectangle dis- 
crepancy of (0, m, 2)-nets lifted to the sphere via the transformation $, which we state in 
the following result. 

Theorem 1. Let P/v — {^o? • • • > x b m -i} Q (0, l) 2 be a (0, m, 2)-net in base b. Let x n = 
(x ni i,x n>2 ) for < n < b m . Let z n = Q{x n ) for < n < b m and Zn = {zo, . . . ,z^-\}. 
Then 

b 2 m 1/ 1\ 1/ 

D^Z N ; S», tl) < ^ ^ + ^ (9 + - j + ^ - 1 

and 

2 b 2 m 1 /9 1\ 1/6 1 

D^{Z N ;S i n)<—— + — ^- + -J+^^-- l 

Proof. The result follows from Lemma [I] and [T41 Theorem 1 and Remark 2] . □ 



46 + 3 
(6+T) 2 



46 + 3 
4(6+l) : 
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Lemma [T] and the result of Roth [35] imply that 

Llog 2 iVj +3 



D N {Z N ;S 2 ,n)>D* N {Z N ;S 2 ,n*) > 



2 8 N 



for all point sets Z^ = {z , . . . , Zjv-i} C S> 2 (where log 2 denotes the logarithm in base 2). 
Thus the construction of Theorem [T] is optimal with respect to the discrepancy defined in 
Definition [3] 

Theorem [T] is not surprising since the transformation $ is area preserving for all rect- 
angles, which follows from ([5]). 



2.2 Uniform distribution on the sphere 

A sequence of iV-point systems {Zn} n>2 on the unit sphere S 2 in M. 3 is said to be asymptot- 
ically uniformly distributed on S> 2 if the exact integral 1(f) = L 2 f da 2 of any continuous 
function / on S 2 integrated with respect to the normalized surface area measure a 2 on 
S 2 can be approximated arbitrarily close (as N becomes large) by means of the equally 
weighted quadrature rules Qn having these Z N as node sets; that is 

Km Q N (f) = 1(f) for every / G C(§ 2 ). 

N— >oo 

This limit relation also states that the weak-* limit (which is defined in such a way) 
of the discrete measure placing a point mass 1/N at every point in Z^ is given by the 
natural measure on § 2 . In other words, the limit distribution (as iV — >■ oo) of the iV-point 
configurations Z N is given by a 2 . Let \Z N \ denote the cardinality of the (finite) set Z N . 
Therefore, an equivalent characterization is that 

lim \Z N HA\/N = a 2 (A) 

N— s>oo 

for every 02-measurable set AC§ 2 (whose boundary has 2-dimensional Hausdorff measure 
0). Informally speaking, any such set A gets a fair share of points as iV becomes large. In 
fact, it is sufficient to consider spherical caps on S 2 . Let C = {C(z,t) : 2 G § 2 , — 1 <t< 1} 
denote the set of all spherical caps. Thus, a natural measure to quantify 'equidistribution' 
of iV-point systems on S 2 is the spherical cap discrepancy 

D(Z N ;S 2 ,C)-= sup 

ccs 2 

where the supremum is extended over all spherical caps. It is well-known that the sequence 
{Z N } N > 2 is asymptotically uniformly distributed on S 2 if and only if D(Z N ;E> 2 ,C) — > 
as iV — > 00. One can show that the following assertions are equivalent (cf., for example, 

my- 

1. The sequence {Z N } N > 2 is asymptotically uniformly distributed on S 2 . 

2. lim^ Q N (f) = 1(f) for every / G C(S 2 ). 

3. limTv^oo <5jv(y«,m) — for all spherical harmonics Y^ m of degree £ > 1 from a (real) 
L 2 (S 2 , o"2)-orthonormal basis {Y^ m }. (This is WeyVs criterion on the sphere.) 



z N nc\ 

N 



<T2(C) 
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4. Hindoo D{Z N - § 2 , C) = 0. 



The spherical cap discrepancy can be estimated in terms of Weyl sums by means of 
Erdds-Turdn type inequalities [] (Grabner [19], also cf. Li and Vaaler [25] ) 



D(Z N 



;S 2 ,. 



^=1 



C2 C 3 

£ L + 



1 £ 



m=l 



1 iV_1 

77 2J Y t,m(Xj) 



N 



3=0 



where L is any positive integer and the positive constants c±, C2, and C3 are independent 
of N, or LeVeque type inequalities ([29], generalizing LeVeques result for the unit circle 

!) 



c 4 



00 Z{d,l) 
£=1 m=l 



iV-1 



i=0 



1/2 



< S 2 ,C) <c 5 



00 Z(d,l) 
£=1 m=l 



JV-1 



A r 



j=0 



l/(rf+2) 



where af.=T(£ — l/2)/T(£ + d+ 1/2) x l/£ d+1 =:bg, for some positive constants C4 and C5 
which are independent of N. The integers 



Z(d,0) 



Z(d, 



(2£ + d-l) 



r(e + d- 1) 
r(d)r(€ + 1) 



denote the number of linearly independent spherical harmonics Y^ m of degree i. 

Instead of using spherical caps as test sets, one can define discrepancy with respect to 
spherical rectangles (as done in this paper) or, as proposed by Sjogren [37], with respect 
to the family of so-called ^-regular test sets. A ^-measurable set A C §> 2 is defined to be 
K -regular if the ovmeasure of the ^-neighborhood (S sufficiently small) of its boundary, 

\z E § 2 : dist(A, z) < 5, dist(§ 2 \A,z)<8\, dist(A, z):= inf \\z - z'\\ , 

is linearly bounded by K 8. Clearly, spherical caps (rectangles) are .fT-regular for some 
K > 0. In pQ Andrievskii, Blatt and Gotz related the discrepancy of a measure with the 
error in integration of polynomials in the following sense. 

Proposition 1. There exists a universal constant Cq such that for every probability mea- 
sure v supported on § 2 and every K-regular set ACS 2 there holds 

K 



\u(A) - a 2 (A)\ < Co inf i — + C(v,n) 
neN n 



where 



C{v, n):= sup 



pdv — j pdo~2 



p polynomial on S 2 , 
degp < qn, \p\ < 1 on § 2 



and q = 30 



^uch inequalities are a generalization of Erdos and Turans result for the unit circle [TBI 117) . 
2 Note that in [1] everything is done in K. . 
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In particular, if v is the counting measure J2xez N induced by the node set Zn of 
an equally weighted quadrature rule Q N on §> 2 , then we have that 



v{A) - a 2 (A) 



z N r\A\ 



I pdu- I pda 2 = — P( z ) ~ j Pd<x 2 



in Proposition [T} 

Let D(Z N ; § 2 , Tk) denote the discrepancy with respect to the family of ^-regular test 
sets and D(Zn]E> 2 ,Q) be the spherical rectangle discrepancy. 

Theorem 2. Let {Zn} be a sequence of N -point configurations with N = b m , m > 1, 
defined by (0,m,2)-nets Pn in base b lifted to the sphere S> 2 by means of Zn = $(P/v). 
Then the following holds: 



(i) linijv^oo D{Zm'-, § , Tk) = for every fixed K > 0. 

(ii) lim Af ^ 00 L'(Z A r;§ 2 , T K ) = 0. 
(Hi) limAf^oo D(Zn, S 2 , Q) = 0. 

Consequently, the sequence {Zn} is asymptotically uniformly distributed on § 2 . 

Proof. Let Pn be an (0, m, 2)-net with N = b m points (m > 1) lifted to the sphere S 2 
and z k = &(x k ) (0 < k < N). Then, using 'cylindrical' coordinates z = (y/l — t 2 z*,t), 
where z* e S 1 and -1 < t < 1, and the decomposition da 2 {z) = (1/2) dt da x {z*) (|25j). 
for every polynomial p on S 2 (satisfying \p(z)\ < 1 on S 2 ) we get 



l N-i . r 1 N-i x 

Qn(p) ~ Hp) = jjJ2p(z k ) ~ pda 2 = J^zZ p ^ Zk ^~2 P( z ) dt 

k —Q J§ 2 "'S 1 u—n J —1 



k=0 



d<7i(«*). 



For each fixed z* G S 1 we define a new iV-point configuration Zn by aligning all points in 
Z N such that they have the 'same z*\ that is z k = (y/l - t\z\t k ) (0 < k < N). Clearly, 
Z N depends on z*. Z N has N points, since the underlying system P N is a (0, ra, 2)-net. 
Indicating the dependence on z* by z k (z*) we write Qn(p) — I(p) = AiQo) + A 2 Q»), where 



i JV ~ 1 



k=0 



N 



k=0 



p(z k (z*)) da x {z* 



r r i N ~ l 1 r 1 

^S 1 [ fc=o J ~ l 



da^z*). 



First, we consider Ai(p). By definition of (0, m,2)-nets every elementary interval 

[a/b d , {a + 1) /b d ) x [a'/6 m - d , (a' + 1) /6 m - d ) , < a < b d , < a' < b m ~ d , 

contains exactly one point of Pn- By construction of Zn, the 2-sphere can be partitioned 
into M x \/N polar zones A\, . . . , of equal sizes I/a/ZV, each containing precisely 
iVg x V^/V points of Z N (N± + • • ■ + N M = N). Picking M heights T\, . . . , r M such that 
the circle Ce with height Ti lies in A^ we move the points in Ai on the circle Ci giving 
them the same height T£ without changing the 'longitudes', that is z' k :—(^/l — rjz* kl Tf) 
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for Zk = (yl — t\z* k , tk) € Ag. The error introduced in this way can be made arbitrarily 
small by increasing N, since the polynomial p is bounded and uniformly continuous on 
§ 2 and the widths of the polar zones uniformly decrease with increasing N. The integral 
J gl p(zk(z*)) d(Ti(z*) averages the polynomial p over the circle on § 2 with height tk- In 
the zone Ai we may use the approximation J gl p((yF 



t}z 



M 



Ai(p) = E 



iV, 



£=1 



1 



2 * 
TtZ 



r T t ))d<Ti(z 

T i ))da 1 (z* 



Therefore 



#i(p), 



E p ( z 'fc) " 

where the error R\(p) goes to zero as N — > oo. Observe, that the square-bracketed 
expression is the error of integration of an equally weighted quadrature rule on the circle 
Ce with integration points z' k for Zk € Ag induced by a horizontal strip of the underlaying 
(0, m, 2)-net Pn for the polynomial p restricted to Ce- Since the extreme discrepancy of 
any (0, ra, 2)-net tends to zero as m — > oo, it follows that the limit distribution of the 
integration nodes is given by <j\ . Therefore the square-bracketed expression tends to zero 
uniformly for all 1 < I < M as iV — > oo. 

Next, we consider A 2 (p). We define the function f(z*; t):=p((\/l — t 2 z*,t)) for \t\ < 1 
and each fixed z* G S 1 , which is uniformly continuous and bounded in t and z*. Thus 



, r, N-i „i 



da^z*). 



The square-bracketed expression is the error of integration for an equally weighted quadra- 
ture rule with integration nodes given by tk = 1 — 2xk, where the XkS are the 2nd- 
coordinates in the (0,m,2)-net Pn, integrating a function from the class {f(z*; •) : z* € 
S 1 }. We can apply Koksma's inequality and use that \f(z*; t)\ < 1 to see that A 2 (p) — > 
asiV-> oo. Thus, for every polynomial p on § 2 (satisfying \p(z)\ < 1 on § 2 ) we get that 
\A(p)\ < |Ai(p)| + |A 2 (p)| -> as iV -> oo. 

Now, we can apply Proposition [TJ Let g = 3. Set fN = J2 z ez N Then to every 
£ > we can chose a smallest integer n such that K/n < e/2 and find an N n such that 
(7(2, u N , n) < e/2 for N > N n yielding 



z N r\A\ 



N 



a 2 (A) 



<Coi^ + C{2,v N ,n)}<C £ 



which holds for every K- regular test set A (K is fixed). Item (i) follows. 

Since spherical caps (rectangles) are X-regular (if'-regular) for some K, K' > 0, the 
Items (ii) and (iii) follow. By the list of equivalent characterizations of uniform distribu- 
tion on S d , the sequence {Zn} is asymptotically uniformly distributed on S 2 . □ 

Concerning the convergence rate of the spherical cap discrepancy, it was an observation 
of Beck [1] that to any iV-point set Zn on S d there exists a spherical cap C such that 

\z N n (7| 



Cl 



AT-3/4 



< 



N 



CT 2 (C) 



and (using probability arguments) there exist iV-point sets Z N on S 2 such that 

\z N nc\ 



N 



(72(C) 



< c 2 N- 3/A ^hgN 



for any spherical cap C. (The numbers q > are constants independent of N.) 
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Remark 2. Thus, the correct order of decay (up to a y/logN factor) of the spherical 
cap discrepancy of a sequence of low discrepancy N -point configurations on S> 2 is given by 
iV- 3 / 4 as N -> oo. 

This is in contrast to the convergence rates (essentially [up to a logarithmic factor] 
1/N) given in Theorem fl] for the discrepancy with respect to spherical rectangles. 



3 Numerical integration on S 2 

We consider now numerical integration of functions over S> 2 in some reproducing kernel 
Hilbert space defined on § 2 . Let P k , k 6 No, denote the Legendre polynomials. We define 
a reproducing kernel 

oo 

K(z,z') = J2XM(z,z')), 

e=o 

for some real numbers Xi > 0. The rate of decay of A^ is related to the smoothness of the 
functions in the associated reproducing kernel Hilbert space. In particular, if A^ x £~ 2s , 
then the reproducing kernel Hilbert space "H s consists of functions with smoothness s, see 
[TT] for more details. 

In [TTI Theorem 5.1] a choice of \i was introduced for which A^ x £~ 3 and for which 
the following reproducing kernel can be written explicitly as 

oo 

K(z,z') = J2X e P e ((z,z')) = 21- \\z-z'\\, 

e=o 

where X is the distance integral 



1=[ [ \\z-z'\\da 2 (z)da 2 (z'). 



This integral can be computed and has value X = 4/3 (see [H]). The coefficients A^ are 
given by 

A = X, A, = X ~[~ 1/2) ^ far/>l, 

where (a)t = a(a + 1) • • • (a + £ — 1) = T(a + £) /T(a) is the Pochhammer symbol. See [11] 
for more details. 

The function K(z, z') is a reproducing kernel of a Hilbert space H 3 ^ 2 over § 2 . In [TTJ it 
is shown that the squared worst-case error for numerical integration in this "H 3//2 is given 
by 

JV-l 

e 2 (Q N ,H 3 / 2 )=l-—J2 \\*h-z t \\. (6) 

k,e=o 

For quadrature rules Qm with node sets maximizing the sum of distances J2k7=o \\ Zk ~ Zi II 



(that is, minimizing the right-hand side in (J6J) which in turn means that the nodes have 
minimal /^-discrepancy), it is known that there are constants c, C > such that 

cAHV 2 < #(Q Ni nW) < CN- 3 / 2 (7) 

for sufficiently large, see again (TTJ Corollary 5.2]. Note that the lower estimate always 
holds. 
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3.1 A lower bound on the worst-case error 

Let Q* N be a quadrature rule whose integration points z*,..., z* N maximize the sum of 
distances Y^kJ=o 11^ — ^11 or > equivalently, have minimal L 2 -discrepancy. (Such points 
always exists by the continuity of the distance function and compactness of the sphere.) 
Then, using the lower bound in ([7]), we have 

1 N-l N-l 

e\Q N ,U i,2 )=T-j f2 £ \\zM-ZiW >X-±- 2 \\zt-zl\\^e\Q* N ,n^)>cN-^ 

k,£=0 k,i=0 

for some positive constant c not depending on N and N sufficiently large. Thus, one 
obtains that c 2 (Qn,'H 3 ^ 2 ) > cN~ 3 / 2 for sufficiently large N for equal weight quadrature 
rule Qn induced by any Appoint configuration. This is in agreement with the results 
obtained in [22] . 

3.2 An upper bound on the worst-case error 

By Stolarsky's invariance principle ^ and representation (|6| we have 

N-l 

e\Q N , U 3 ' 2 ) = X - — J2 \\z k ~ zi\\ = 4 [D 2 (Z N ; § 2 , C)] 2 . (8) 

k,e=o 

Since the spherical cap discrepancy provides an upper bound for the L 2 -discrepancy via 
D 2 {Z N ;E> 2 ,C) < V2D(Z N ;S 2 ,C), one has necessarily e 2 (Q N ,n 3/2 ) — )• as A^ — )• oo for 
a sequence of asymptotically uniformly distributed Appoint configurations on S> 2 . This 
applies in particular to digital nets on the sphere of the type considered in Theorem [2j 
A natural question is if such nets yield the same order of convergence as the worst-case 
error over the unit ball in % 3 ^ 2 over S 2 provided with the reproducing kernel K(x, y) = 
(8/3) — 1 1 x — y || has (which is the same as of the optimal L 2 -discrepancy on § 2 (see ffify 



as suggested by the numerical results in Section 3.3 



Theorem 3. Let {Z^} be a sequence of N -point configurations on S> 2 defined by point 
sets P/v Q [0, l] 2 lifted to the sphere S 2 by means of Zn = $(Pjv)- Then the equal weight 
quadrature rule Qn associated with Zn satisfies 

e 2 (Q N ,n 3/2 ) < (^= + 2^) D*(P N ;[0,l) 2 ,n*). 

Proof. Notice that f g2 \\z — z'\\ da 2 (z f ) does not depend on z G § 2 . Thus, defining the 
function 

1 N ~ X 

A(w):= / \\w — z\\ da 2 (z) — — \\w — z k \\ , w e § 

J& k=0 

one obtains 

N-l „ „ - N-l 



(9) 



i £ A(z e ) = [ f \\z - z'\\da 2 (z)da 2 (z') -±^J^\\ Zk - Zi \\ = e 2 {Q Nl U 3 / 2 ) 

£=0 J® 2 kl=0 
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and therefore 



e 2 (Q N ,U 3/2 ) < max A(z e ) < sup A(z). 



Note that the sum in (J9J) can be seen as the potential (with respect to the distance kernel) 
of the discrete measure associated with Zn at w. Thus, A(zi) gives the deviation from 
the leading term of the asymptotic expansion (as iV — > oo) of this potential at Zg. 
Let 

g(a,r) = \\w - $(a,r)||. 
Since, by (|5]), the transformation $ is area- preserving, it follows that 

/ \\w — z\\ d(T2(z) = / / g(a, r) da dr. 
Js 2 Jo Jo 

Thus, for any fixed w e § 2 , we can view A(w) as the integration error of integrating the 
function g using the points Pn- By the Koksma-Hlawka inequality, this integration error 
is bounded by the star-discrepancy D*(P N ; [0, 1] 2 ,1Z*) of Pjy times the total variation of 
the function g in the sense of Hardy and Krause. 

The variation of g in the sense of Vitali is given by 

V^ 2 \g) = su P J2\ 5 (9,J)l 

v 

where the supremum is taken over all partitions of [0, l) 2 into subintervals J = [a%, 02) x 
[61,62) an d S(g,J) = g(ai,bx) — g(ai,b 2 ) — #(02, &i) + g(a>2, 62)- The variation of the 
one-dimensional projections is given by 

M 

V {l \g) = sup ^2\g(x k ,l)-g(xk-i,l)\ 

0=x <x±<—<x M =l 

and 

M 

V^ 2 \g)= sup ^2\g{l,y k ) - gO-,Vk-i)\. 

0=J/Q<J/1< — <J/Af=l k= l 

The total variation of g in the sense of Hardy and Krause is then given by 

V(g) = V (1 \g) + V^(g) + V^ 2 \g). 

Let w = (0, 0, ±1). Then 5(g, J) = for all rectangles J and therefore V^ 1 ' 2 )^) = 0. 
Further, $(a, 1) = (0,0,-1) for all < a < 1 and therefore V {1) (g) = 0. Finally, 
^ {g) = 0) — (7(1,1)| = 2 because of the monotonicity of g(l,-)- Thus we have 
V(g) = 2. 

Let w 7^ (0,0, ±1). Then there are u, v G (0, l) 2 such that w = $>(u,v). Again, we 
have V^(g) = 0. Further, we have V^ 2 \g) < 2y2, where we have equality if u = 1 and 
v = 1/2. 

Now we consider V^ 1 * 2 ^^). The variation V lyl,2 \g) does not change by changing u 
because of the rotational symmetry of g about the polar axis. We may therefore assume 
without loss of generality that u = 1/2. First, let v = 1/2. The mixed derivative 

a2 x2 (1 — 2t) ( 1 + a/(1 — r) r cos(27ra) ) sin(27ra) 
dadr { ' ' drda { ' ' r- , rr. r~ A 3 ' 2 



y/(l-r)T (2 + -t)t cos(27ra)V 
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is finite for < a < 1 and < r < 1 with $(a,r) ^ w. In particular, its sign does 
not change for (a, r) in the interior of one of the four quadrants Ti = [0, 1/2] x [0, 1/2], 
% = [0, 1/2] x [1/2,1], T 3 = [1/2, 1] x [0, 1/2], f 4 = [1/2, 1] x [1/2, 1] of the square. For a 
subinterval J with J contained in the interior of one of these quadrants, one can use 

%■/) = J^g^(x)dx = VOL(J)-^-(xj) for some x j G J (10) 

to determine the sign of S(g; J). For a subinterval J which shares an upper or lower 
boundary with the square, the quantity S(g;J) reduces to a difference of two function 
values, since g remains constant for r = 0, 1. In this case one can use monotonicity of 
the function g along horizontals (decreasing towards a = 1/2). For J with w as one 
corner, one can still deduce if either 5(g; J) < or 5(g; J) > 0. It suffices to consider 
partitions V w defined by horizontal and vertical lines including the lines with t = 1/2 and 
a = 1/2. From our observations on the sign of S(g; J) we conclude that the expression 
J2jeP w \$(9> J)\ can ^ e reduced to two summands for each quadrant (and by symmetry) 

V (1 ' 2 \g) = £ \S(g, J) | = -2(0(0, At') - g(l/2, At') + - 2) + 2(g(0, At') - g(l/2, At')) 
Jev w 

= +2(2 - + (7(1/2, 1 - At") - g(0, 1 - At")) +2(^(0, 1 - At") - g(l/2, 1 - At")), 

where At' and At" denote the "heights" of the first and last row of subintervals of V w . 
Cancelling terms, we arrive at 

V^ 2 \g) = 8 for («,„) = (1/2,1/2). 

Therefore, for the case v = 1/2, we can use the Koksma-Hlawka inequality to obtain 

\A(w)\ < D* N (P N )V(g) < D* N (P N ) (8 + 2y/2) . 



Now assume that < v < 1/2 (the case 1/2 < v < 1 follows because of symmetry). 
Here, the mixed derivative is given by 

(j c)'~ (j 

'a, t) = (a, t) = — 16t (1 — v) v (1 — 2t) i?(cos(27ra)) sin(27ro;) [g(a, t)] _ , 



dadT dTda 
where the linear form R is given by 

v - t + (l - 2v) (l - t)t 



R(x) = x + 



v / (1-v)v{1-t)t (1 - 2t) ' 

Hence the mixed derivative is again finite for < a < 1 and < r < 1 with $(«, t) ^ w. 
Furthermore, it vanishes on the lines a = 0, 1 and on the line a = 1/2 (except when 
t = 1/2, where it is undefined). It does not vanish on the line r = 1/2. Finally, it 
vanishes on the set (cf. Figure [TJ 

T:= {(a, t) : (1 - 2r) i?(cos(27ra)) = 0,0<a<l,0<T<l}. 

(Note that the mixed derivative is singular everywhere along the lower and upper bound- 
ary of the unit square which are mapped to the poles of the sphere.) The set T is symmetric 
in the sense that (a, t) £ Y if and only if (1 — a, t) G V . From the observation 

, , I 0+(l-2v)(l-v)v 

R(x)\ =x+—^ /— = x + l 

v ' w=v 1 -v)v{l-2v) 
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it follows that T and the line r = v intersect only in the point (1/2, v). We observe that 
the following three relations are equivalent: 

\R(x)-x\ < 1, 
(v -t)(v- 3vt + 3vt 2 - r 3 ) < 0, 

r > v and r 3 + 3vr < 3vr 2 + v. 

From the last pair we obtain the equivalence 

\R(x)-x\<l if and only if v < r < v + (1 - v) 2/3 v 1/3 - (1 - v) 1/3 v 2/3 =:r v . 

We conclude that T is the graph of a curve contained in the strip v < r < t v which is a 
function r(a) (assuming the value t v at a = 0, 1 and v at a = 1/2) and is a function a(r) 
when restricted to either the left or right half of the unit square. We record that t v < 1/2 
ifO < v < 1/2 and t v = 1/2 if v = 1/2. Moreover, the vertical line test also shows 
that the function r(a) is strictly monotonically decreasing towards a = 1/2. It is also 
continuous and continuously differentiable as can be seen from an implicit differentiation 
of -R(cos(27ra:)) = (the factor (1 — 2r) can not be zero in our setting): 

T f (a) = — An sin(27ra)— ^— — — — ^ . - — , where r = r(a). 

v ; v 'v-3vt + 3vt 2 -t^ + 3{t-v){1-t)t y ' 

(The denominator is a strictly monotonically increasing function on the interval (v, 1/2).) 
Using the relation i?(cos(27ra;)) = along T in order to substitute cos(27ra) in g(a, r), we 
obtain g restricted to T and its first derivative in the forms 

t(o) — v , , v 2 T'(a) 



g r {a):=g{a,r(a)) = 2\j T ^_ , g' T {a 



1 - 2vg v {a) 

Finally, we observe that g is decreasing along horizontal lines as a tends to 1/2, since 

^-(a, t) = -8n^(l-v)v(l-r)r Sm{2n ^ , (a, r) G [0, l] 2 , (a, r) ± (1/2, v) 
ool g{a,r) 

and g is decreasing along vertical lines in the strip f<r<l/2asr tends to v, since 

2t)x. 



dg_ 
dr 



(a, r) = 2 y ,J , H(x) = 1 - 2v + J y - J — (1 

g{a,T) y (1 - r)r 



Similar as in the case v = 1/2, the curve T and the vertical line a = 1/2 divide the unit 
square into four regions Ti, T2, T3, T4. In the interior of each region the mixed derivative 



has the same sign. By (10), the sign of 8 (g; J) is the same for all subintervals contained 
in a fixed region T k . Suppose V w is a partition of the unit square determined by a grid 
which is symmetric with respect to a = 1/2 and also contains the vertical a = u(— 1/2) 
and the horizontals r = v and r = r v . (It suffices to consider this type of partitions.) The 
horizontals r — v and r = t v subdivide the regions further into a rectangle T' k without 
the curve V and a part T£ with Y as part of the boundary. By symmetry (cf. Figure [T]) 

V (1 ' 2 \g)= \5(9,J)\ = 2(g(0,Ar')-g(l/2,Ar'))-2^ 5(g,J)-2 £ S(g, J) 

jnr=0 
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u u 

Figure 1: Contour plot of g(a, r) with curve T and scheme to compute V^ 1,2 \g). 
+ 2 l%> J )l+ 2 E S(g,J) + 2j^ 5(g,J) + 2(g(0,l- At") - g(l/2, At")) , 

J£P W j_ JSP™, JeVw, 

JCT1UT4, JCTf, JCTl 

Jnry0 jnr=0 

where At' and Ar" denote the "heights" of the first and last row of subintervals of V w 
and in the sums one excludes the J's from the first and last row. Using cancellation 
at "interior" nodes, we obtain 

V^\g) = 2 fo(0, Ar') - ^(1/2, Ar')) + 2 (g(0, v) - - g(0, At') + 0(1/2, Ar')) 

-2^%,J) + 2 l%^)l + 2^%,J)+2(9(0,l-Ar")-9(l/2,Ar")) 

Je"P™, J^P WJL _ JeVw, 

JCT{', JCT1UT4, JCT^', 

jnr=0 Jnryo jnr=0 

+ 2 (0(0, t v ) - 0(1/2, r„) + 0(1/2, 1 - Ar") - 0(0, 1 - Ar")) 

= 2g(0,v) -2 J) + 2 ]T |%J)| + 2 ^ %J) + 2( S (0,r„)-#,r t )). 

Je"P™, JeP^ JsP™, 

JCTf, JCTiUT 4 , JCT^', 

jnr=0 Jnry0 jnr=0 

Let v = y m < y m +i < ■ ■ ■ < y m +K+\ = t v denote the r-values of the grid defining V w in 
the strip v < r < t v . To every y m +k with < k < K — 1 there exists a maximal x^ m+fc of 
the a-values defining the grid such that the subinterval [x(_ k _ x , y m +k) x l x e m+k , Vm+k+i) 
does not intersect T. For consistency reason we set l m +K = 0, so that Xi m+K = 0. Then 

K-l 

- E 5 ^9,J) = E W * 2W*+i) ~ P(^ m+ fc'^+jfc+i) +fi , (^ m+fe '?/™+fc) - P(°>2/m+fc)) 
JePu,, fc=o 

JCTf 7 , 

jnr=0 

K-l 

= g(0,y m +i<) -g(0,y m ) - E (9(x£ m+k ,y m +k+i) - g{xi m+k ,y m+k )) ■ 

k=0 
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Since g is decreasing along horizontals in the strip t><r<l/2asr goes to v, one gets 

- ^ J "> - 9(0,y m + K ) - g(0,y m ) < g(0,r v ) - g(0,v). 

JCTf, 

jnr=0 

Similarly, to every y m+k with 1 < k < K there exists a minimal x Pm+k such that the 
subinterval [x Pm+k , y m +k) x [x Pm+k+1 , y m +k+i) does not intersect T. For consistency reason 
we set p m to that index with x Pm = 1/2. Then 

K 

6 (9,J) = Yl (9( x Pm+k^ym+k) - g(l/2,y m+k ) + g(l/2,y m+k+1 ) - g(x Pm+k , y m+k +i)) 
Je-P™, k=\ 

jnr=0 

K 

= g{l/2,y m+K +i) -g{l/2,y m+1 ) - ^ (g(x Pm+k ,y m+k+1 ) - g(x Pm+k , y m+k )) 

k=l 

< <7(l/2, Vm+K+i) - g(i/z,y m +i) < g(l/2,y m +K+i) = g(l/2,r v ). 
As an intermediate result we have 

V^ 2 \g)<2g(0,v) + 2(g(0,T v )-g(0,v)) + 2 £ |% J)| + 2<?(l/2, r v ) 

J&T WJ _ 
JCT1UT4, 

Jnry0 

+ 2((/(0,r„)- 5 (l/2,r t ,))=4y(0,r t ,) + 2 ]T 

JCT1UT4, 
Jnry0 

Using triangle inequality and monotonicity along horizontals we get 

K Pm+k-l 

Y l% J )l=E Y \9( x j^V'm+k) - g(x j+1 ,y m+k ) + g(x j+1 ,y m+k+1 ) - g(x j ,y m+k+1 )\ 

J&P W , k=0 j=£ m+k 

JCT1UT4, 
Jnry0 

K Pm+fc-1 K Pm+fc-1 

<J] ^ {g(xj, y m +k) - g(x j+ i,y m +k)) ^ (g(xj,y m+k+1 ) - g(x j+1 , y m+k+1 )) 

fc=0 j=£ m+k k=0 j=£ m+k 

K K 

= Y {g( x e m +^ym+k) - g(x Pm+kJ y m+k )) + ^ (g(x tm+k ,y m+k+1 ) - g(x Pm+kJ y m+k+1 )) 

k=0 k=0 
K 

^ 2 J2 (9rk + J - gr( x Pm+ J) ■ 

k=0 

In the last step we used that the left side of a horizontal segment of successive rectangles 
of the covering of T is below and the right side is above the curve T by construction, cf. 
Figure [lj Recall, that gr is strictly monotonically decreasing as a tends to 1/2. Since 
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9r(ze m +k-i) > 9r(x Pm+k ) for < k < K - 1 by monotonicity of g T , we have 

K 

Y \5(g,J)\<2{ (gr(x e J - g r (x p J) + Y M^wJ ~ flrOw*-i)) 

Jep w ^_ fc=i 
JCT1UT4, 
Jnry0 

+ Y {9r( x em+k-i) ~ 9r{x Pm+k )) } = 2^(0) - ffr(l/2) + ^ (sr(^ m+ *-i) - ffrO^+J) }• 
fc=i fc=i 

By construction , „ . > x v ,„> xp , „ „ > x D , „ , > • • • > Xp , , > x v = 1/2. So, 
because of monotonicity of g-p, we increase the right-hand side of the estimate above when 
including terms terms for the gaps [x Pm+k ,X£ m+k _ 2 ]. Thus 

Y \S(9,J)\ <4(^(0)-sr(l/2)) = 4sr(0), 

JCT1UT4, 
Jnrv0 

where the expression in parentheses is the total variation of (monotone) one-dimensional 
function gr restricted to [0, 1/2]. Putting everything together, we obtain 



V {l > 2 \g) <4<7(0,r„) + 2 Y \8(9,J)\<49(0,T v ) + 8g T (0) = 12g T (0) = 24J^ 



, . 2v 

JCT1UT4, 
Jnry0 



It can be shown that the right-most side above as function of v is strictly monotonically 
increasing on [0, 1/2) with 

,. ,'r v -v 1 
lim 



V-+1/2- V 1 - 2f v^' 
We can use the Koksma-Hlawka inequality to obtain 

|A(w)| <D*(P N ;[0,l] 2 ,TV)V{g) < D*(P N ; [0, l] 2 , TV) ^ + 2^ for < v < 1/2. 

This concludes the proof. □ 

Digital nets [15], [31] are explicit constructions of point sets in [0, l) 2 with D*(Pn; [0, l] 2 , TV) 
0{N~ X log AT). By mapping them to the unit sphere S 2 one also obtains explicit construc- 
tions of points on S> 2 with favorable discrepancy. We have the following corollary. 

Corollary 1. Let b > 2 and m > 1 be integers and let N = b m . Let {Z^} be a sequence 
of N -point configurations on S 2 defined by digital (0,m,2)-nets Pn € [0, l] 2 lifted to the 
sphere S 2 by means of Z N = $(P N ). Then the equal weight quadrature rules Qn associated 
with Z N satisfy 

e 2 (Q N ,n 3/2 ) = OiN' 1 log N) as 00. 
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3.3 Numerical results 



We consider now the worst case error ([6]) for quadrature rules defined by digital nets 
based on a Sobol' sequence, see [39J. For efficient implementations of Sobol' sequences see 
[21 El [23] . Figure [2] shows the squared worst-case error e 2 (Qn, 7~L 3 ^ 2 )- The numerical results 
suggest that the squared worst case error 6 2 (Qn, % 3//2 ) converges with order 0(N~ 3 ^ 2 ) as 
N ->• 00. 

Note that the first point of the Sobol' sequence (or any digital (0, 2)-sequence for that 
matter) is always (0,0). Since $(0,0) = (1,0,0), this point gets mapped to the North 
Pole. On the other hand, no point of the digital sequence gets mapped to the South 
Pole. This might not be a desirable feature. To remedy this situation one can randomize 
the (0, 2)-sequence using a scrambling algorithm, see [271 El]- In this case the sequence 
in [0, l] 2 is still a (0, 2)-sequence, but the point (0,0) only occurs with probability 0. 
Numerical investigation of scrambled Sobol' sequences yield similar results as the one 
shown in Figure [2] 




Figure 2: The dashed lines show N~ 3 ^ 2 and (9/4)iV _3//2 , and the curve shows the squared 
worst-case error e 2 (Qjv, % 3//2 ), where the quadrature points are a digital net mapped to 
the sphere. 



The numerical results lead us to the following conjecture. 

Conjecture 1. Let x ,x 1 ,... e [0, l) 2 be a (0,2) -sequence and let z ,z 1 ,... G S 2 be 

the corresponding points on the sphere obtained by using the mapping $. Let QnIJ) = 
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m 


N = 2 m 


e 2 (Q N ,H 3 / 2 ) 


AT" 3 / 2 


N^e\Q N ,U 3 l 2 ) 


1 


2 


6.2622e-01 


3.5355e-01 


1.7712 


2 


4 


2.1149e-01 


1.2500e-01 


1.6920 


3 


8 


8.1448e-02 


4.4194e-02 


1.8430 


4 


16 


3.5091e-02 


1.5625e-02 


2.2459 


5 


32 


8.0526e-03 


5.5242e-03 


1.4577 


6 


64 


2.6309e-03 


1.9531e-03 


1.3470 


7 


128 


9.4336e-04 


6.9053e-04 


1.3661 


8 


256 


3.4501e-04 


2.4414e-04 


1.4132 


9 


512 


1.3374e-04 


8.6316e-05 


1.5495 


10 


1024 


4.6029e-05 


3.0517e-05 


1.5083 


11 


2048 


1.8846e-05 


1.0789e-05 


1.7468 


12 


4096 


6.4670e-06 


3.8146e-06 


1.6953 


13 


8192 


1.7873e-06 


1.3486e-06 


1.3252 


14 


16384 


5.6815e-07 


4.7683e-07 


1.1915 


15 


32768 


1.9912e-07 


1.6858e-07 


1.1811 


16 


65536 


6.3194e-08 


5.9604e-08 


1.0602 


17 


131072 


2.4122e-08 


2.1073e-08 


1.1447 


18 


262144 


9.1906e-09 


7.4505e-09 


1.2335 


19 


524288 


3.7001e-09 


2.6341e-09 


1.4047 


20 


1048576 


1.3068e-09 


9.3132e-10 


1.4032 



Table 1: Numerical results: The worst-case error obtained when using digital nets over 
Z 2 lifted to the sphere S 2 . 

~k T,n=o f( z n)- Then we have 

e 2 (Q N , n 3/2 ) = 0(N~ 3/2 ) as N ->■ oo. 

In other words, a (0, 2)-sequence lifted to the 2-sphere via the mapping <£> achieves the 
optimal rate of convergence of the worst-case integration error in U 3/2 . By Stolarsky's 
invariance principle, the conjecture also implies that a (0, 2)-sequence lifted to the 2-sphere 
via the mapping <£> achieves the optimal rate of decay of the spherical cap L 2 -discrepancy. 

Acknowledgement: The first author is grateful to the School of Mathematics and 
Statistics at UNSW for their support. 
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